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Abstract 

The formalism of two coupled Dirac equations within constraint instant form dynamics is used 
to study the nucleon-nucleon (NN) interaction. The salient features and the final Schrodinger type 
equation is given. Explicitly energy dependent coupled channel potentials, for use in partial wave 
^) ' Schrodinger like equations, with nonlinear and complicated derivative terms, result. We developed 

£N| , the necessary numerics and study np and pp scattering phase shifts for energies 0-3 GeV and the 

^q' deuteron bound state. The interactions are inspired by meson exchange of tt, T), p, u> and a mesons 

, for which we adjust coupling constants. This yields, in the first instant, high quality fits to the 

' Arndt phase shifts 0-300 MeV. Second, the potentials show a universal, independent from angular 

momentum, core potential which is generated from the relativistic meson exchange dynamics. Ex- 
trapolations towards higher energies, up to TLab = 3 GeV, allow to separate a QCD dominated short 
I , range zone as well as inelastic nucleon excitation mechanism contributing to meson production. A 

' local and/or nonlocal optical model, in addition to the meson exchange Dirac potential, produces 

agreement between theoretical and data phase shifts. Third, the 1 5o, 3 Po and 3 Pi partial waves 
elicit a fusion/scission, for Ti, a b < 1 GeV, and a fusion/fission, for T^ab > 1 GeV, mechanism for 
, intermediate dibaryon formation. 

o : 

O ■ 1 Introduction 

The formalism of coupled two-body Dirac equations, within constraint instant form dynamics, is used 
'-^ , to study the nucleon-nucleon (NN) interaction. This particular approach for two spin 1/2 particles 

^ \ was developed by Crater, Van Alstine, Long and Liu [21 03 [U El • They define a Poincare invariant 
interaction in terms of scalar, pseudo scalar, vector etc. interactions with the implication that they 

J> . satisfy certain compatibility conditions [7j. This approach yields in its final form explicitly energy 
dependent coupled channel potentials for use in partial wave Schrodinger like equations We followed 

5_i ■ and re-derived their expressions up to a certain point and developed our own numerics to study np 
and pp scattering phase shifts for < T^ ab < 3 GeV. The comparison with recent data makes use of 
SM00 and SP03 GWU/VPI SAID phase shift solutions 0. 

The NN interaction is described within the paradigm of exchange mechanism involving w, p, to a 
and other mesons exchanges to make up what we call the Dirac potential. A comparison with 
most recent experimental data, GWU/VPI SP03 phase shifts, requires the adjustment of coupling 
constants and a regularization of the short range interaction domain. For energies above pion produc- 
tion threshold 280 < Ti, ab < 3000 MeV we added a phenomenological complex optical model potential 
to the Dirac potential. This addition brings the theoretical S-matrix in perfect agreement with the 
experimental data S-matrix from GWU/VPI and, more importantly, permits the identification of some 
predominant reaction paths. 

1 Presented at XXII International Workshop on Nuclear Theory, Rila Mountains, June 16-22 2003, 
E-mail: geramb@uni-hamburg.de 
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The coupled two-body Dirac equations, combined with the meson exchange model, yield the ap- 
pearance of a repulsive, practically hard core, potential independent of partial wave. The universal 
core radius has a value r c = 0.5 ± 0.025 fm. This core radius is independent of a nucleon substructure. 
It depends only on masses, in particular of the exchanged mesons, and the full relativistic treatment 
of the NN system. This feature is not present with equal distinctness in any of the current NN best fit 
potentials of np and pp data [HI ITT]. For purpose of comparison, we show results of the Argonne 
AV18 potential [TO]. 

The fitting process of coupling constants uses data in the sub-meson-production domain, < 
Tiab < 280 MeV, of np and pp partial wave phase shifts. For > 280 MeV, single and double 

intrinsic nucleon excitations, A(3, 3) and other low excited hadrons, as well as simple and complex 
reactive meson productions contribute. This is well known and demands beyond NN a more complex 
coupled channels problem to solve. We curtail the problem to NN scattering using an optical model 
potential (OMP) addition ^3- Despite of complicated inelasticities, selection rules of angular 
momentum, isospin selection and the complex energy dependences some of the partial waves show that 
the real phase shifts 5{T), are very well reproduced (extrapolated) by the Dirac potential alone. Most 
clearly, this is realized in the 1 S'o, 3 -Po an d 3 -Pi channels and < 1100 MeV. We have realized this fact 
before |12l I13j but wish now to support more convincingly the case of an intermediate fusion/scission, 
for T^ai, < 1 GeV, and a fusion/fission, for T^f, > 1 GeV, mechanism in which two nucleons change 
briefly into a compact dibaryon with subsequent decay back into two nucleons and mesons. 

In Fig. 1 we show an intuitive and guiding scheme which distinguishes interaction domains as 
function of separation between the two nucleons. This scheme is in accordance with coupled channels. 
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Figure 1: NN scattering and reaction scheme for < 3 GeV. 
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2 Theoretical background 



Relativistic quantum mechanics demands a Poincare invariant formulation. This implies unitary repre- 
sentations of the Poincare group as transformations for state vectors. One distinguishes ten generators 
for translations and rotations 

P M (4) and «V(6) = -«V (1) 

The associated commutation relations yield the Poincare algebra. Subsets of these generators are 
associated with subgroups, which correspond to 3D hypersurfaces in Minkowski space. The kinematic 
subgroups are instant form x° = 0, light front form x° + x 3 = and point form x ■ x = a 2 > 
with x° > jTj. Within the fundamental variables, one distinguishes simple kinematic variables and 
complex Hamiltonians. 

Crater, Van Alstine and their collaborators El El El 121 treat the nucleon-nucleon problem with 
two coupled Dirac equations. Each of the two free nucleons satisfies 

(ij^df, - m)i/f = and (p 2 + m 2 )^ = 0. (2) 

They relate to a representation in terms of Todorov's variables 

x = x\ — X2 relative distance 
P = pi + P2 total momentum 

ui = \/—P 2 invariant mass, in CM. total energy 



(P' z + w 2 )* = 0, P = in CM. 
x\_ = {rf^ + P u P^)x^ transverse coordinate 



P» 

UJ 



9 9 9 9 9 9 

UJ + ml — 771-2 00 + m 2 ~ m l 

61 = Yuj ' 62 = 

777,0, = — - — - relativistic reduced mass 



9 9° 

uj — m| — m<2 



UJ 

relativistic relative energy on reduced mass 



2uj 

e 2 pi - eip 2 eiP e 2 P 
p= , P-p = 0, pi = \-p, p 2 = p 

UJ UJ UJ 

2 2 2 u 4 + (m 2 -m^) 2 -2uj 2 (m 2 +ml) 

k =p = b (uj, mi, 7772J ~ 



4w 2 

2 2 2 

L - m 1 = e 2 ■ 

In particular ^-matrices are related to 7-matrices 



222222 

e u, ~ m u = e l - m l = e 2 - m 2- 



^ = ^ = diag(-l,l,l,l), ^ = 7^757^, 5 = M/^75, 

ei = {rf» -p u p^/p 2 )e^ 

Several anticommutator relations hold, viz. 

[e»,<n+ = -rT, [9 5 ,e 5 \ + = -i, 
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The free equations are 

(p- 6 + m6 5 )?p = and (p 2 + m 2 )^ = 0. (3) 
The generators of the Lorentz group J^ u contain angular momentum and spin 

Jfiv — Lfjtv + S^i/ 

L ^v = 2 (%P» + P»% - Q»P» - P/J-Qv) 
% % 

Each nucleon moves in the field of the other nucleon. Ultimately, the interaction is described by a 
meson exchange model. 

For two particles with spin 

Sio = PiOi + mi8 5i , i = l,2 

the commutator 

[<5iOi<52o] = 

vanishes strongly. 

2.1 Example with scalar interaction 

A scalar interaction changes the mass into a mass operator Mi = rrii + S. The Dirac equations are 

Slip = (+p ■ 6i + eiP • 9i + Mi0 5 i)V> = 

s 2 ip = (-p ■ e 2 + e 2 P ■ e 2 + M 2 e 52 )ip = o 

[S 1 ,S 2 ]-ip = -i (dMi ■ 01052 + dM 2 ■ 9 2 9 51 )ij + 0. 
The commutator does not vanish in general but vanishes through a third law condition 

d{Ml (x± ) - Ml ( x± ) ) = and M\ (x± ) - M 2 2 (x x ) = m\ - m\ . (4) 
Using a hyperbolic parameterization 

M\ = 777-1 COsh L + 777,2 sinh L 

M 2 = 777-2 cosh L + mi sinh L 

gives, in terms of a single invariant function L = L(x±), a compatible representation of the scalar 
interaction. The two body Dirac equations, with a scalar interaction, are 



Siif) = (+0i • p + ei0i • P + Mi 05i - idL ■ 2 05205iJ i> = 

S 2 4> = (-02 • P + e 2 2 • P + M 2 52 + idL ■ O^O^ i> = 

dM x dM 2 



where dL 



M 2 Mi 
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The Dirac constraint operators satisfy 
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(Sf - S$) ij> = - - (pi + ml - p\ ~ ml) i> = -P ■ p i> = 0. 



(5) 



Crater, Van Alstine and collaborators generalized their hyperbolic representation to apply for any 
interaction being Poincare invariant, viz. 



Scalar: 



A, 



L #51 #52 



L 



-0\ with Oi = -751752 



Timedike Vector: 



Spacedike Vector: 



Pseudo-scalar: 



Aj = JP ■ 9 1 P ■ 9 2 = 2 - = /3i/3 2 ^Oi 



A g = g e 1± -e 2± = o 3 ^ = 7i±72±|ci 



c c 

Ac = £\— = -151I52O1 — 



and four others of which we shall not make any use. With such sum of interactions 

A = Aj + A L + A g + A c , 

the coupled system of equations 

<SiV = (+G6 l ■ p + Erf! ■ P + Mi 51 + £ (6> 2 • dG0 3 + J0 2 - LO x + C£{) ) 



(6) 



S 2 il> 



-G6 2 ■ p + E 2 6 2 ■ P + M 2 52 - £ {Oi ■ d0O 3 + J0 2 - LO x + C£ x ) 



are fully defined by masses, energies and single particle quantities 

Mi = mi cosh(L) + m 2 sinh(L) 
M 2 = m 2 cosh(L) + mi sinh(L) 

E\ = ei cosh( J) + e 2 sinh( J) 
E 2 = e 2 cosh( J) + e\ sinh( J) 

and G = e G . 

Finally an elaborate Pauli reduction yields coupled Schrodinger type equations. <r, p, u, it and 77 
exchanges specify the interactions 



L{x±),J(x ± ),C(x±),G(x±). 



(7) 
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We use 
and 



V = E X M 2 + E 2 M 1 



B 2 = E 2 - M 2 = E 1 - M? 



j x — ivi 1 

The final stationary Schrodinger type equation 



P 2 ~i 



2gl _E 2 M 2 + Ml E 1{j + Ly 



V 



i(J-L)' 



(((71 • f)(<T 2 • P) + (<T2 • f)((Ti • p)) 



l -v 2 g-\g' 2 - \{C + J - L)'(-C + J - L)' + 



- & (J + L) 



+ (o"l • cr 2 ) 



2 2 



1 E 2 M 2 + M\Ei n i i -f t r N / l„ w ia' 1(-C + J-L)' 



L • (cJi + cr 2 ) 



-g'u + l)' - -g'C - 

7 2 2 r 2 



' 1 E 2 M 2 + M\E\ 

"2 p ( J + L ) 



+ 



L • (<Ti — a 2 ) 1 £ 2 M 2 - ^1-^1 
r 2 P 

L ■ (o"i x a 2 ) i M 2 E X - M X E 2 
r 2 V 



(J + L)' 
(J + L)> 



+ (cJl • f)(fJ 2 • f) 



■\v 2 {-c + j - l) - \v 2 g - g'i-c + g' 2 + ^-g' 

11 2r 



= e- 2 ^B 2 \6 + ) 



(8) 
(9) 



(10) 



can be treated with well known techniques, in particular partial wave expansion, to find angular 
momentum (with spin, isospin and angular momentum) dependent radial Schrodinger equations. 



3 Dirac Potentials and Partial Wave Phase Shifts 

Partial wave expansion separates spin, angular and radial parts to yield coupled or uncoupled radial 
Schrodinger type equations 

-<f>"(r, k) + Vi(r)</>'(r, k) + V o (r)0(r, k) = k 2 (j)(r, k). (11) 

The appearance of a first derivative term requires solving a system of first order equations which 
is numerically not favorable. However, the Numerov algorithm is favorable and popular for radial 
Schrodinger equations but applies only to second order, without first derivative, equations. A factor- 
ization of the solution (j)(r, k) by the ansatz 

4>(r, k) := g(r, k)f(r, k) 
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yields first order equations 

g'(r,k) = —Vi(r)g(r,k) with lim g(r,k) = 1 (12) 

or second order equations 

9"(,r, k) = (\v{(r) + - A V?(r)) g(r, k) with lim g(r, k) = 1. (13) 

The other factor satifies the second order equation 

-/"(r, fc) + F(r)/(r, fc) = A; 2 /(r, fc), lim f(r, k) ~ ^(-/T(r, fc, r?) + fc+(r, fc, r?) S(fc)). (14) 

Ricatti Hankel or Coulomb functions h (r,k,rf) and the S-matrix S(k) determine the asymptotic 
solutions. We identify the Dirac potential with 

V D (r) = g'Hr, k) Qt?(r) - ^V{(r) + F (r)J g(r, k) (15) 
and additionally centrifugal and Coulomb potentials 

V FC (r) = g~\r, k) (^+3. + Bj^j ^ fc) . (16) 
We use the convention of E] 

P = P\+ P2 total momentum, 



\/—P 2 invariant mass, in CM. total energy, P = in CM., 

w 2 + m 2 — m\ <^ 2 + m 2 , — 
61 = 2^ ' 62 = 2^ ' 

?n w = — - — - relativistic reduced mass, 
u 

u> — m\ — m<2 



relativistic relative energy, 

t2Pi ~ eiP2 e x P , e 2 P 
P= , P\ = hp, P2 = P, 

UJ UJ OJ 



2 2 uj a + (m\ - ml) 2 - 2o> 2 (m 2 + mf) 2 2 
k =p = — 2 = e w -m u = e x - m l = e 2 - m 2 , 

2 



6 6 

r](k) = — 7—dpp Coulomb parameter 

fv 



This is to be compared with standard non-relativistic NN potentials JO]- The recent work by Liu and 
Crater [2j elaborates on other methods to determine phase shifts from Eq. 1)11(1 . 

The model specification for L(x±), J(x±), C(x±) and G(x±) follows Liu and Crater model I [5] to 
specify scalar 

q * e ~ m ° r ( \ 2 e ~ ma ° r 2 e- m *> r 
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pseudo scalar 

C = (n • r 2 )^- + ^- ^- . 

lo r uj r uj r 

and vector 

p-m p r p-m^T p-m^r 



A = (n • r 2 ) 9 ; + s£ + , 



interactions. All Yukawa form factors are regularized with a normalized Gaussian 



For S < we use 



and for 5 > 



— 'VTL X 

dx 3 - e -(^) 2 /« 2 with a = 0.14142 fm. (17) 

x 



smh(L) = 1 + , 18 



Ml = m\ + G 2 (2m w S + S 2 ) 

Mf = m 2 + (-2(2^^ + 

<9Mi <9M 2 



M 2 Mi 



This models tt, rj, p, u, 5 and a exchanges. In Figs.|3J0J[2]and|l)]are shown Dirac potentials for three 
values, T^ab = [0.1, 1,2] GeV (red lines). In comparison are shown the results of the popular Argonne 
AV18 potential (blue line) ^Uj. The remarkable feature of the Dirac potentials is their universal 
repulsive core with r c ~ 0.5 fm. The only exception is the 3 PF± channel where the ansatz of repulsion 
turns, surprisingly, into a short range attraction. 

In Figs. [7| and |H1 are shown the phase shifts of SM00 (green), SP03 (blue) and theoretical results 
(real Dirac potential solutions (red), real Dirac potentials with complex OMP added are coinciding 
with the data of SP03 blue lines). 

The Dirac instant form dynamic yields partial wave spin, isospin and energy (a channel) dependent 
NN potentials V^{r : T) to which we add a local or nonlocal optical model potential whose strengths 
are fitted to data 

V(r) = V°(r, T) + [U a (T)g(r) + iW a (T)g{r)) (19) 

or 

V(r,r') = V°(r,T)6(r - r') + \<j){r) > [U a (T) + iW a (T)] < 0(r')|. (20) 

Since local/nonlocal potentials imply similar results, we restrict ourselves here to the local optical 
potential and reference the more general and nonlocal case |12j . 

A remark about phase shift convention appears necessary. Consider the partial wave radial equation 

-/4'(r, k) + V a (r)f a (r, k) = k 2 f a (r, k) (21) 

with stationary 

/° ~ j(kr) + n(kr)K(k) (22) 
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or physical boundary conditions 



/+ ~ h~h~(kr) + h+(kr) S(k)). (23) 
1% 



S- and if-matrices are related to phase shifts |14j . There exist several conventions to represent the 
S-matrix in terms of phase shifts. The Arndt and Roper Jl] (GWU/VPI) convention uses S- and 
K-matrices 

S= (1 + iK)(l - IK)' 1 , with K = i(l - 5)(1 + S)~ l = ReK + ilmK. (24) 

Re if corresponds to a unitary S-matrix and phase shifts 5 ± (k) and e{k) defined by 

q(^K\-( cos(2e)exp(2i(5-) i sm(2e) exp(i{5~ + S+))\ 

o{nejx) y isin (2e)exp(i(8- +5+)) cos(2e) exp(2i<T ) J" 1 j 

Absorption phase shifts, p^ and /i, are related to 

/ tan 2 p~ tan tan p + cos u\ ,„„.. 

ImA = _ , ^2+ • I 26 ) 

ytan p tan /r cos /x tan /r y 

Single channels simplify to if = tan 5 + i tan 2 p. 



4 The Nucleon-Nucleon Optical Potential 

The notion of an optical model is useful in cases when the S-matrix is not unitary and flux disappears 
into open inelastic or reaction channels. The optical model is often expressed in terms of a complex 
and energy dependent potential where the imaginary part effectively describes the loss of flux without 
specification of the inelastic channels. A less popular alternative to a complex optical model potential 
is the introduction of pseudo channels. Here we follow the optical potential approach |12j . 

The source of the BB channel (dibaryon formation) is the NN core-domain whose transition is 
mediated by a delta-function or a narrow Gaussian function. Intrinsic nucleon excitations are also 
mediated by a narrow Gaussian in the core domain. Inelasticities are either generated by coupling 
BB (dibaryon) to asymptotic many body final states, composed of two nucleons and mesons, or decay 
of XY, composed of one or two intrinsic nucleon excitations, into asymptotic many body final states. 
Within the inner core region r < r c NN and XY wave function components vanish. The meson exchange 
Dirac potential, which is described by the NN Dirac instant form dynamics, should ultimately be limited 
to r > r c ~ 0.5 fm in its effect. This constraint eliminates the need for regularization of the short range 
Dirac potential and boundary conditions are automatically generated by the 5(r — r c ) NN <-> BB 
transition potentials. This proposal is demonstrated in Fig. |2 The strengths and location of 5-function 
interactions are boundary conditions which are to be determined by BB and XY models. Herein lies 
the essential point of our method. Dirac potentials play only the role of a shield which prevents us from 
seeing the naked refinement surface of hadronic QCD dynamics - it recalls the P-matrix formalism. A 
realization of the full coupled channels problem is in progress. 
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REGULARIZATED NOT REGULARIZATED 




R [fm] R [fm] 



Figure 2: Dirac potential (red), left figure regularized and right figure not regularized, for the np 1 So 
channel, T^ab 0.1, 1 and 2 GeV, showing a long range OPEP tail, an attractive pocket ~ 0.75 fm and a 
core repulsion with r c = 0.5 ± 0.025 fm (blue). The regularized and not-regularized potentials show, in 
all channels, only small differences. Regularization does not affect the conclusion drawn about the core 
geometry generated but helps to keep numbers reasonable near the origins. Also inserted are Gaussian 
form factors g{r — r$) ~ exp(— (r — ro) 2 /a 2 ) which are used with the optical model (thin line curve). 
In this figure tq = r c = 0.5 fm and a = 0.2 fm. 



Without specification of details, the coupling scheme has the following structure 

7~Lbb + ^BB^( r ~ r c) + Vbb confined dibaryon, NN meson decay (27) 

Vj?&*(r-r c ) + V$%g(r - r c ) + H NN + NN elastic (28) 

g(r - r c ) + H X Y + Vxy NN meson decay. (29) 

Below pion production threshold, Ti a b ~ 280 MeV, the NN S-matrix is unitary for all practical 
purposes. Above this energy excitation of A(3,3) resonances is the predominant mechanism. It is 
obviously present in the NN 1 D2, 3 F^ and 3 Pi ? 2 channels. Isospin conservation suppresses a coupling 
to NA in the np T = channels. NN scattering, for energies below 3 GeV in general, show a comparable 
to nucleon-nucleus scattering weak and smoothly energy dependent coupling to inelastic channels. A 
perturbative treatment of inelastic and reaction channels with DWBA methods is thus strongly favored. 

A key issue for all secondary applications of NN scattering is a high quality reproduction of the 
elastic NN scattering channel. Inverse scattering methods are useful for this purpose. These methods 
use the experimental data in form of partial wave phase shifts as input and determine the optical model 
potential as a correction to a theoretically defined and numerically realized reference potential. 
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Figure 3: Dirac potentials (red) for singlet 5 = 0, 
channels. Differences between np and pp potentials 
for pp. The AV18 channel potential is shown as blue 
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Figure 7: Single channel, pp [0,3] GeV, rep [0,1.2] GeV SM00 (green) and SP03 (blue) real phase shift 
data 5(T) and theoretical curve based soly upon Dirac potential (red). The full potential, Dirac 
and OMP, describes the data SP03. Of particular interest are 1 Sq, s Pq, 3 P\ channels in which the 
real phase shift of the Dirac potentials coincide with data for [0,1.1] GeV but diverge above 1.1 GeV. 
Obvious are effects due to A(3, 3) coupling in other channels. SM00 (continuous solution (green), single 
energy solutions, open green circles, with error bars). SP03 (continuous solution (blue), single energy 
solutions, full blue circles, with error bars). Theoretical results (real Dirac potential solutions, full red 
line, real Dirac potentials with complex OMP added are coinciding with the data of SP03 blue line). 
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Figure 8: Coupled channels pp [0,3] GeV, np [0,1.2] GeV SM00 and SP03 real phase shift data 5(T) 
and theoretical curve based soly upon Dirac potential (red line). The full potential, Dirac and OMP, 
describes the data SP03. SM00 (continuous solution, green, single energy solutions, open green circles, 
with error bars) and SP03 (continuous solutio, blue, single energy solutions, full blue circles, with error 
bars) and theoretical results (real Dirac potential solutions, red line, real Dirac potentials with complex 
OMP added are coinciding with the data of SP03 blue line). 
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Figure 9: Available absorption phase shifts p(T). Dirac potentials generate no absorption and OMP 
parameters Uq(T) and Wq(T) are adjusted to reproduce the continuous energy solutions of p(T). 
SM00 (continuous solution, green, single energy solutions, open green circles, with error bars) and SP03 
(continuous solution, blue, single energy solutions, full blue circles, with error bars) and theoretical 
results (real Dirac potential solutions, red, real Dirac potentials yields no absorption but with a complex 
OMP added the results are coinciding with the data of SP03 blue line). 
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4.1 1 S , 3 ^o and 3 P 1 Channels 

NN phase shift data, see Figs.[2|and|ni show in almost all channels and for Ti ab > 280 MeV a compli- 
cated energy dependence and deviations from the Dirac potential predictions. Exceptional cases are 
the 1 5o, 3 -Po an d 3 -Pi channels, contained in Fig.[7[ which show a practical perfect reproduction 

Soata = <Wac, for < T Lab < 1100 MeV (30) 

by the Dirac potential alone. However, the absorption 

PData^O, whereas p Dirac = for 280 < T Lab < 1100 MeV. (31) 

This demands an optical potential U(r) + iW(r) that leaves the real phase shifts 5(k) unchanged but 
generates an absorption p(T) > for 280 < T Lab < 1100 MeV. 
The optical potential solutions are complex 

~fr + (V° + U)f r - k 2 f r = +Wfi 
~f- +{V° + U)f i -k 2 f l = -Wf r 

and the boundary conditions, S = 5(V D , U, W) and p = p(V D , U, W), are transcendental functions of 
potentials and OMP adjustable parameters. 

A sensible solution is U(r) = and W{r) = W5{r — tq) as optical model interaction in Eq. (|19j) . 
We used a delta-function and/or a narrow normalized Gaussian, see Fig.|2 

w(r) = {w G N(r ,a )exp(-(r - r ) 2 /«§) f , 

\W 5 5(r-r ) 

with U = 0, starting at the origin with 

/r(0)=0, f r (h) =h i - l+l \ fi(r) = for < r < r . (33) 

The following conclusions are drawn: The phase-shifts <5(T), p(T) for 280 < T^ ab < 1100 MeV imply an 
optical potential at the surface of the repulsive core for 1 5o, 3 -Po and 3 P\ partial waves. Intermediate 
dibaryons are practically not formed, the BB channel is realized by a dibaryon fusion/scission picture as 
shown in Fig. 1111 The BB dibaryon quark dynamic is reduced to energy dependent complex boundary 
conditions at the core radius permitting meson production. The meson exchange mechanism is not 
valid inside the core radius. Caveat, at this stage of our work, we integrated from the origin through 
the core region realizing a small real wave function at the core radius. This, in connection with the 
Gaussian OMP form factor, causes a real and imaginary part OMP Ug(T) ^ and Wg(T) ^ 0. 
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Figure 11: Left, caused by the Pauli exclusion principle for a six-quark dibaryon and TL a b,NN < 
1100 MeV suggests a futile ringing of the nucleons and a suppression of dibaryon formation. It gives 
the impression of a fusion/scission mechanism. Right, the formation of dibaryons, at sufficient high 
energy, is governed by medium to longer ranged quark-gluon flux tubes with fusion into a six-quark 
hadron sized dibaryon with sequential decay. This inspires a fusion/fission mechanism. 

In Fig.EH are shown the np OMP (normalized Gaussian, ro = r c = 0.5 fm and ao = 0.2 fm) 
strengths values. The crucial center of baryon NN and BB transition radius is ro = 0.5 ± 0.025 fm. 
At higher energies > 1.1 GeV, the transition surface becomes more and more faded, washed- 

out and translucent when the energy of dibaryon states matches the total energy of the NN system. 
Intermediate short lived dibaryons J BB = + , 0~, l 1 are formed, see Fig, [Til 

We estimate, from the phase behavior in these three channels, the total energy (lowest mass) of a 
dibaryon system m BB = 2400 ± 150 MeV and a width T > 150 MeV. 

Coupling to dibaryons is realized for Ti a b > 1100 MeV and the fusion/scission picture may change 
gradually into a fusion/fission picture as shown in Fig. 1111 
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